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We  report  experimental  data  which  trace  the  time  development  of  electric  discharge  channels  in  air 
and  which  demonstrate  the  turbulent  cooling  of  such  channels.  These  data  provide  qualitative  confir¬ 
mation  of  the  model  proposed  and  used  by  Hill,  Rinker.  and  Wilson  to  calculate  the  production  of 
nitrogen  oxides  by  lightning.  We  outline  an  analytical  treatment  which  identifies  assymmetries  in  the 
pressure  and  density  gradients  or  the  discharge  channel  as  a  significant  source  of  vorticity .  The 
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r_> equations  from  which  the  vorticity  strength  and  mixing  time  scale  may  be  calculated.  We  briefly 
describe  detailed  simulations  with  which  we  have  calibrated  the  theory.  Finally,  we  combine  the 
experimental  data  with  our  calibrated  formulas  to  estimate  the  convective  mixing  rate  in  the  case  of 
lightning.  We  obtain  a  rate  oft300cm3/sec  per  cnr^air  in  the  return  stroke  channel  after  pressure 
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CONVECTIVE  COOLING  OF  LIGHTNING  CHANNELS 
1.  Introduction 

Hill  et_  al.  (1980)  have  recently  proposed  a  "channel  heating"  model 
to  predict  the  production  of  atmospheric  nitrogen  oxides  (NO^)  by  lightning. 
In  this  model,  the  return  stroke  deposits  energy  over  ~100us,  heating  the 
lightning  channel  to  ~30000K  within  a  few  microseconds.  As  energy  is 
deposited,  the  heated  region  (i.e.,  the  channel)  expands  and  compresses  the 
atmosphere  surrounding  the  channel  to  produce  a  shock  wave.  Subsequently, 
the  shock  wave  decouples  from  the  channel  as  the  return  stroke  current  drops 
to  small  values.  The  channel  continues  to  expand  until  pressure  equilibrium 
is  reached  with  the  surrounding  air  at  a  temperature  of  ~3000K.  The  rate 
at  which  the  air  inside  the  channel  cools  is  important  to  the  calculation  of 
NO^  production  by  lightning.  Hill  et_  al.  (1980)  calculated  the  production 
rate  based  on  the  assumption  of  "turbulent"  or  convective  cooling,  in  which 
the  ambient  air  is  mixed  with  the  air  inside  the  channel  at  a  prescribed 
constant  rate.  The  authors  used  a  detailed  air  chemistry  computer  code  to 
simulate  cooling  at  several  rates,  the  lower  limit  of  which  was  based  on  the 
rate  of  cooling  by  thermal  conduction  in  the  absence  of  convective  mixing. 

They  left  open  questions  regarding  the  cause  of  the  convective  cooling,  the 
efficiency  of  mixing,  and  the  length  of  the  time  interval  between  pressure 
equilibration  of  the  channel  and  the  onset  of  mixing. 

In  this  note,  we  report  experimental  data  which  directly  confirm  quali¬ 
tative  aspects  of  the  above  model  and  which  demonstrate  the  turbulent  cooling 
of  discharge  channels.  We  then  outline  a  preliminary  analytical  treatment 
which  identifies  asymmetries  in  the  pressure  and  density  gradients  of  the 
discharge  channel  as  a  significant  source  of  vorticity,  which  in  turn  causes 

the  ambient  air  to  mix  into  the  channel.  Our  theoretical  analysis  results 
Manuscript  submitted  January  12.  1081. 
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in  equations  from  which  the  vorticity  strength  and  the  mixing  time  scale  may 
be  calculated.  To  calibrate  this  analytic  model  and  test  its  predictions, 
we  have  performed  detailed,  two-dimensional  numerical  simulations  which  will 
be  described  briefly  and  compared  with  our  analytic  results.  Finally,  we 
combine  the  experimental  data  with  our  calibrated  formulas  to  estimate  the 
convective  mixing  rate  in  the  case  of  lightning. 

2.  Experimental  Studies 

We  have  recently  performed  experimental  studies  of  the  propagation  of 
electric  discharges  through  air  and  the  dynamics  of  the  resulting  channels 
(Greig,  et  al. ,  1978).  In  Fig.  1,  we  show  Schlieren  photographs  of  discharges 
which  were  15  cm  in  length.  The  discharges  were  produced  by  a  Marx 
generator  with  an  erected  voltage  of  ~250  kV.  The  current  had  a  peak  value 
of  ~10  kA,  a  period  of  ~3,i,s,  and  a  duration  of  ~6y,s.  While  this  is  a  much 
shorter  pulse  than  that  occurring  in  a  lightning  return  stroke,  we  expect 
the  qualitative  behavior  to  be  similar  after  the  shock  decouples  from  the 
channel  because  air  response  in  rotational  flow  is  relatively  slow.  We 
have,  in  fact,  performed  preliminary  experiments  using  parameters  similar 
to  those  of  a  lightning  return  stroke.  We  could  not,  however,  obtain  a 
complete  sequence  of  photographs  similar  to  Fig.  1  because  the  channel 
size  exceeded  the  capability  of  our  optical  system. 

In  Fig.  1,  we  see  that  a  hot,  smooth,  curved  channel  with  a  radius  of 

~1.4  cm  has  formed  within  8p,s  of  the  initiation  of  the  discharge.  The 

accompanying  shock  wave  has  already  decoupled  frcji  the  channel  and  appears 

as  a  sharp  line  at  the  boundary.  By  30us,  the  shock  is  easily  identifiable 

and  propagates  out  of  th®  field  of  view  at  a  time  somewhat  greater  than  lOQ^s 

after  discharge  initiation.  At  approximately  10C>4,s,  the  interior  temperature 

18  -3 

of  the  channel  is  ~5000K,  the  gas  density  is  ~30  cm  and  the  electron  density  is 
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~  LO  cm  .  (These  values  are  the  measured  conditions  for  a  comparable  discharge 
in  a  pre- ionized  channel.)  Th-‘  heated  channel  remains  roughly  stable  up  to  that 
point.  By  250j.-s,  however,  instability  has  begun  to  distort  the  channel,  as 
evidenced  by  density  fluctuations  at  the  channel  boundary.  The  photograph  at 
900  shows  that  the  distortions  become  more  pronounced  as  cooler  air  at  the  edges 
of  the  channel  is  pulled  toward  the  axis.  The  average  radius  has  also  increased 
to  ~  1cm,  indicating  an  increase  of  ~  100X  in  the  volume  of  the  channel  due  to 
entrainment  of  the  surrounding  air.  As  the  Lime  from  discharge  initiation  increases, 
smaller  scale  (turbulent)  structure  appears  and  by  10ms,  the  channel  has  begun  to 
disappear.  We  estimate  the  final  channel  radius  to  be  ~6-7cm,  indicating  a  volume 
increase  of  ~  20  times  the  original  volume.  Studies  of  similar  discharges,  which 
have  been  guided  along  a  path  designated  by  lnser/aerosol  interaction,  have  produced 
nearly  identical  results.  Channel  cooling  and  the  onset  of  turbulence  were  some¬ 
what  more  rapid  for  such  "laser-guided  discharges . "  This  result  follows  directly 
from  the  analytic  theory  described  in  the  next  section. 

S.  Theoretical  and  numerical  studies 

2 

Our  theoretical  treatment1'  attributes  the  "turbulent"  (convective)  cooling 
phenomenon  to  asymmetries  between  tire  gradients  of  the  existing  pressure  and 
density  distributions.  These  asymmetries  will  generate  vorticity  as  the  channel 


expands,  according  to  the  equation 

,  5 
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a  significant  residual  vorticity  exists.  This  vorticity  is  responsible  for 
mixing  of  ambient  air  with  the  hot  channel  gas.  We  should  point  out  that  other 
mechanisms  for  generating  residual  vorticity  in  the  heated  channel  might  exist. 
For  example,  rapid  movement  of  the  discharge  current  axis  may  occur  as  a  result 
of  magnetic  forces  present  when  the  current  is  nonnegligible.  Such  movement 
could  cause  sufficient  displacement  of  the  surrounding  air  to  produce  some  long¬ 
term  mixing  motion.  We  will  reserve  consideration  of  such  phenomena  for  a  future 
more  detailed  paper. 

We  have  identified  several  types  of  asymmetry  which  might  be  relevant 
to  lightning  and  the  above  discharge  experiments: 

(1)  Two-dimensional  distortions  of  the  return  stroke  channel  from  a 
circular  cross  section, 

(2)  Two-dimensional  asymmetries  from  return  stroke  displacement  off 
the  axis  of  the  leader  stroke,  and 

(3)  Three-dimensional  distortions  in  the  return  stroke,  such  as  the 
curvature  of  the  channel  axis  shown  in  Fig.  1. 

If  we  define  a  cylindrical  coordinate  system  with  the  z-axis  along  the 
average  direction  of  the  channel  axis  (for  channel  sections  which  do  not 
undergo  large  changes  in  direction),  we  find  that  all  three  types  of 
asymmetry  will  produce  a  nonzero  value  of  This  will  cause  mixing  in 

the  (r.d)nlanc.  In  addition,  for  three-dimensional  asymmetries ,  will 
be  non-negligible,  producing  mixing  in  the  (r,z)  plane  (a  vortex  sheet). 

Fig.  2  shows  that,  for  sufficiently  short  channel  sections,  all  three  types 
of  asvmmetrv  may  be  represented  by  a  displacement  in  the  (r,~ )  plane, 
which  causes  the  rotational  symmetry  about  the  z-axis  to  be  broken.  In 
our  derivations,  we  will  define  the  x-axis  to  be  parallel  to 
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Fig.  2  —  Schematic  diagrams  of  three  types  of  return  stroke  asymmetry  which  generate  vorticitv: 
< a>  two-dimensional  distortions  (e.g..  elliptical  channel),  (b)  displacement  from  axis  of  leader 
channel,  (c)  three-dimensional  distortions.  The  degree  of  asymmetry  may  be  expressed  in  terms 
of  the  quantity  X0  as  indicated.  In  Fig.  2a  and  2b,  we  show  only  the  portion  above  the  x  —  z 
symmetry  plane. 
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With  this  convention,  rig.  3  shows  the  flow  field  in  the  x-y  plane  for 
cases  b  and  c  in  Fig.  2  after  the  return  stroke  channel  has  reached  pressure 
equilibrium  with  the  ambient  air.  This  flow  pattern  is  approximately  equiv¬ 
alent  to  that  of  a  vortex  filament  pair  having  average  strengths  of  t  »■.  and 
located  respectively  at  (x,  t  y ) .  (The  flow  pattern  for  the  elliptical  channel 
in  Fig.  2b  is  equivalent  to  the  superposition  of  two  vortex  filament  pairs. 

For  simplicity  we  will  discuss  only  a  single  vortex  filament  pair.)  The 
average  strength  a  is  given  by  (Batchelor,  1967) 


and  the  coordinates  are 

X(T) 


.  Uj 


r  (  X  ,  V  .7  ) 

z 


,;2(T) 


I  dy  /  dx  x 
0  * 


x  iz(x,y,T) 


and 


d>:  v  ;  ( x  ,  v ,  t  ) 


( 0/ 


(4) 


(5) 


We  obtain  i  (x,v,t)  by  integrating  (1)  over  tne  time  interval  (0,t),  where 


t  =  0  is  the  time  of  discharge  initiation. 


The  vortex  filament  pair  will  migrate  as  a  unit  at  a  speed  of 


v  .  -  >■.  /4ttv 

pair  z 


(t) 


in  a  direction  parallel  to  the  x-axis.  The  fluid  velocity  along  the 
x-axis  is 


v  (x,0)  = 


:/T(x‘ 


_0 

o 


(7 ) 
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Vx  (x,o) 


Fig.  3  —  Symmetric  distribution  of  localized  vorticitv  appears  as  two  extended  vortex  filaments 
separated  by  a  distance  2y.  Flow  velocity  induced  by  these  vortices  is  shown  along  the  x  axis  half-" 
way  between  the  filaments.  The  mixing  time  rmix  is  defined  here  as  the  time  it  takes  a  fluid  element 
at  -Rj  to  reach  Rj  and  hence  effectively  bisect  the  hot  channel. 
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If  we  define  a  characteristic  mixing  time  scale  t  .  to  be  the  time  interval 

mix 

over  which  a  fluid  element  travels  along  the  x-axis  from  one  "edge"  of  the 
channel  to  the  opposite  "edge,"  we  may  use  (7)  to  obtain 


7mix  "  2*VV'3  +y‘)/l*J-v 


(a) 


where  R  is  the  approximate  radius  of  the  return  stroke  channel  when  pressure 
equilibrium  is  reached.  We  may,  therefore,  characterize  the  mixing  time  scale 
if  we  know  the  vorticity  strength  and  the  average  displacements  of  the  equiv¬ 
alent  vortex  pair  from  the  (x,z)  symmetry  plane. 


We  hive  performed  the  integral  in  (3)  by  using  (1)  in  con j an :  cion  vith 
the  following  assumptions: 

(1)  We  treat  the  return  stroke  as  an  instantaneous  pressure  pulse  of 
finite  size  with  a  given  radial  profile  and  with  a  total  energy  deposition 
equivalent  to  that  of  a  return  stroke.  The  assumption  of  instantaneous  depo¬ 
sition  fits  closely  the  experimental  situation  of  Fig.  1.  At  the  same  time, 
we  do  not  parallel  the  work  of  Hill,  (1971,  1979),  who  permits  energy  to  be 
deposited  over  a  time  interval  ~10CUs.  However,  the  formalism  which  we  have 
developed  is  adequate  for  this  latter  case  as  well. 

(2)  The  time  t  in  (3)  is  equal  to  the  time  required  for  the  return 
stroke  channel  to  achieve  pressure  equilibrium  with  the  surrounding  atmosphere. 
Tor  times  t  >  t  a  residual  vorticity  exists,  and  tnic  vcrticity  is  responsible 
for  convective  cooling  of  the  lightning  channel. 
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(3)  During  expansion  of  the  return  stroke  channel,  tne  fluid  outside 


where  the  boundary  of  the  heated  region  is  R(t)  and  the  velocity  of  the  bound 
is  U(t)  =  R(t).  This  assumption  appears  to  be  reasonable  since  the  shock  wave 
carries  away  only  about  10%  of  the  total  energy  deposited  by  the  return  stroke. 
Further,  this  assumption  is  more  accurate  when  the  energy  deposition  is  slower. 
Notice  that  the  acceleration  is  radial  and  has  values  strictly  appropriate  only 
if  the  effects  of  the  nonaligned  density  gradients  on  the  driving  expansion 
flow  are  small . 

(4)  We  approximate  the  integrand  in  (3)  in  terms  of  its  value  at  time 

t=t  ,  at  which  the  expansion  flux  R  l'  =  R(t  )  U(t  )  is  a  maximum, 
m  m  m  m  m 

Under  the  above  conditions,  we  obtain  for  the  residual  vorticity 


*  VT>  *  VW  *n  <VV  '  <VV  V  V  <10> 

for  the  average  z-component  of  the  vorticity  strength  where  t  ►*.  The 
quantities  R^  and  R^  are,  respectively,  the  initial  average  radius  of  the  neat-- 
(return  stroke)  region  and  the  average  radius  following  expansion,  and  Pa  it  tn 
ambient  density.  The  quantity  is  the  density  at  the  center  of  the  return 
stroke  channel  for  asymmetry  types  1  and  3  (two  and  three  dimensional  dis¬ 
tortions  of  the  return  stroke  channel).  For  asymmetry  type  2  (return  stroke 
displacement  off  the  leader  channel  axis),  kq  is  the  density  at  the  center 
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of  the  leader  channel.  The  dimensionless  form  factor  f  is  a  complicated 
integral  which  contains  geometric  effects,  detailed  hvdrodynamical  inter¬ 
actions  and  information  about  the  leader  and  return  stroke  pressure  and 
density  profiles.  For  asymmetry  type  2,  we  have  evaluated  the  form  factor, 
fos  for  several  different  pressure  and  density  profiles  and  have  found  that 
f  0  is  not  strongly  profile  dependent.  The  values  of  I  f 0  I  vary  smoothly  from 
zero  at  Xq/R^  =  0  to  a  peak  of  <  0.5  at  X^/R^  ~  1.  The  values  then  decrease 
somewhat  more  slowly  to  zero  as  X^/R^  increases  beyond  a  value  of~l.  Similar 
behavior  occurs  for  the  other  asymmetry  types.  For  rough  estimates  of  j ;.  j 
a  value  of  I f !~  ~r  is  permissible,  and  we  mav  substitute  c  for  U  ,  where 
is  the  speed  of  sound  at  the  axis  of  the  return  stroke  when  the  current  flow 
begins  ( t  =  0) . 

We  have  performed  preliminary  detailed  simulations  of  the  problem 
using  the  FAST2D  computer  code  (Boris,  1977)  to  validate  and  calibrate 
the  approximate  analytic  model  developed  above.  In  addition  to  accounting 
for  shocks  properly,  which  the  theoretical  model  does  not  do,  the  simulations 
are  capable  of  describing  the  late-time  motions  and  profiles  as  modified  by 
the  induced  vorticity.  As  above,  we  have  treated  the  return  stroke  as  an 
instantaneous  pulse  of  finite  size  and  equivalent  energy  rather  than  depositing 
the  energy  over  a  characteristic  time  interval  (~10(Ls;.  These  simulations 
have  indicated  that  the  above  analytic  model  and  the  resulting  formula  given 
in  (10)  will,  in  most  cases,  provide  estimates  of  which  differ  from  the 
simulation  results  by  less  than  301.  X'c  have  discovered  that  only  modest 
asymmetry  (Xq~  R^/10)  is  required  to  generate  considerable  mixing  and  that 
for  a  wide  range  of  values  of  X^/R^  (0.1  i.  X^/R^  u  1.5)  the  residual  vorticity 
and  mixing  rate  vary  by  less  than  30%. 


Preliminary  investigations  of  the  asymmetries  of  type  (3),  i.e.,  tnree- 
dimensional  distortions  of  the  return  stroke  profile,  indicate  that  (10;  holds 
for  as  well.  For  cases  in  which  more  than  one  type  of  asymmetry  are  present, 
the  cumulative  effects  should  increase  the  mixing  rate,  and  the  channel  should 
cool  more  quickly.  This  would  explain  the  faster  cooling  rates  for  electric 
discharges  guided  by  preformed  laser  channels,  since  some  displacement  of  the 
discharge  axis  from  the  axis  of  the  preformed  laser  channel  would  be  expected 
as  well  as  three-dimensional  distortion  of  the  electric  discharge  channel. 

Ue  can  now  analvzc  the  mixing  shown  in  Fig.  1,  which  uc  assume  in  this  case 

to  be  due  primarily  to  three-dimensional  curvature  with  a  value  of  X^/R^  ~  0.25. 

We  have  R^  ~  0.1  cm  (measured  from  open  shutter  photographs  of  the  discharges), 

R.  ~  1.4  cm,  P  /Pn  ~  10,  If. I  ~  y  and  U  ~  c  ~  3.5  x  10^  cm/s.  From  (10)  we 
1  “  0  3  3  ms 

-  4  2 

calculate  k  ~3.5  x  10  cm  /s.  Using  (8),  and  y  ~  0.8  R  (from  our  simulations), 

Z  J. 

we  obtain  Tmix  =>-  430  ,j,s.  From  Fig.  1,  we  notice  that  violent  perturbations  begin 
to  occur  sometime  between  ~250p,s  and  ~900p.s.  This  is  consistent  with  the  defi¬ 
nition  and  value  of  t  .  , 

mix 

From  Fig.  1,  we  also  find  that,  for  the  laboratory  experiments,  reasonably 

4 

complete  mixing  occurs  at  a  time  ~  10  ^s,  which  is  ~  ^®Tmix  following 

discharge  initiation.  Referring  to  our  discussion  of  the  laboratory  discharges 

in  Fig.  1,  we  see  that  the  volume  increase  of  the  channel  IV  , (during  the 

cool 

interval  At  ,  )  is  also  ^  20  V.,,  where  V,,  is  the  initial  volume.  We  have 
cool  0  0 

found  that,  in  general  for  convective  cooling,  (Appendix  A), 

"^COOl  .  ^COol  ,  ,,N 


The  ratio  on  the  right  hand  side  of  (11)  depends  only  on  the  temperature  of 
the  channel  before  appreciable  mixing  occurs  and  the  ambient  temperature.  For 


our  laboratory  discharges  and  lightning,  these  temperatures  are  approximately 


the  same.  Thus,  we  assume  that  for  lightning  _t  .  ~  20  .  ,  from  vnich  we  mav 

cooj.  mix 

obtain  rough,  estimates  of  the  time  requirec  for  "complete"  cooling  o:  a  return 

stroke  channel  and  of  the  average  mixing  rate.  For  lightning  we  assume  (Kill 

et  al .,  1980)  R,,  ~  1cm,  R,  ~  16cm,  10 ,  j  f  j  *'•  4  ,  and  l'  ~  c  ~  SxlO'cm ■  s  . 

—  —  U  1  co  i1  J  m  s 

-  A  2 

From  (10)  we  obtain  *•  ~  58x10  cm“/s .  Using  (8)  with  y  ~  0.8  R^.  we  find 

-3  -3 

7  .  ~  3.4x10  s.  Thus  we  expect  "complete  mixing"  to  occur  in  ~  68x10  s  after 
mix 

energy  deposition  by  the  return  stroke.  If,  as  in  Fig.  1,  the  volume  of  the 

channel  increases  to  20  times  its  original  value  during  that  time,  we  will 
3  3 

have  ~  20cm  of  ambient  air  mixed  with  each  1cm  of  air  initially  in  the  channel. 

3 

This  gives  us  a  mixing  rate  (F^  in  the  notation  of  Hill  et_  al . ,  1980;  of  ~  300cm 
3 

per  cm  of  air  originally  in  the  channel,  which  falls  in  the  upper  end  o:  the 
range  considered  by  Hill  et^  al.  Because  we  have  not  treated  the  cumulative  effec 
of  the  various  asymmetry  types  and  because  some  of  the  numbers  which  we  nave  used 

3 

are  only  rough  estimates,  we  feel  that  300cm  /s  should  be  treated  as  an  order  ot 
magnitude  estimate  of  FQ.  Given  our  results,  the  estimates  of  Hill  et  al ■  ( 1 9 S u > 
for  global  NO^  production  by  lightning  appear  to  be  quite  reasonable. 
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APPENDIX  A 


In  this  appendix,  we  give  a  heuristic  proof  of  the  relation  (11) 
for  convective  cooling  of  discharge  channels: 


iV  ,  At 
cool  cool  , 

V  ~  T 

0  mix 


(ID 


We  define  At  ,  to  be  the  time  interval  during  which  the  channel  cools 
convectively  to  ambient  temperature;  ^coo^  is  the  increase  in  channel  volume 
due  to  entrainment  of  ambient  air  over  the  time  interval  At  , ;  V„  is  the 
channel  volume  before  appreciable  mixing;  and  Tmix  has  been  defined  as  the 
time  required  for  a  fluid  element  to  move  from  one  3ide  of  the  channel  to 
the  opposite  side.  We  assume  that  Che  channel  is  cylindrical  with  radius 
and  lengtn  Z  and  that  a  vortex  filament  pair  exists  at  x  *  0,  v  ~  :  R^. 


/ 

/ 


1 


i 
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Fluid  elements  on  different  streamlines  move  at  different  speeds,  and  thus 
adjacent  fluid  elements  on  different  streamlines  will  soon  have  different 
temperatures  and  densities.  Thermal  conduction  and  molecular  diffusion  will 
occur  rapidly  between  adjacent  fluid  elements  since  their  displacement  is 
small,  and  we  therefore  have  enhanced  cooling  over  the  same  situation  (hot 
channel  in  cool  ambient  air)  without  convection.  In  addition,  the  experi¬ 
mental  data  (see  Fig.  1)  indicate  that  turbulent  mixing  occurs  as  the 
surrounding  air  is  pulled  into  the  hot  channel.  These  factors  should  combine 
to  produce  rapid,  efficient  mixing, and  the  continued  entrainment  of  ambient 
air  will  cause  the  channel  to  expand.  To  obtain  the  volume  increase  over 
the  interval  ^‘tcoo^«  we  may  therefore  compute  the  volume  of  ambient  air 
flowing  into  the  channel  per  unit  time.  This  is  given  by  (Batchelor,  pp. 
75-76) 


dV  f 

dt  “  /  £  *  §&  (1A) 

SURFACE 

ABZ 

where  v  is  the  fluid  velocity  and  the  integral  is  taken  over  the  surface 
defined  by  the  channel  axis  Z  and  the  dashed  line  AB  in  the  figure.  Given 
the  definition  of  T  ,  ,  we  may  approximate  v  by: 


mix 


The  area  of  the  surface  ABZ  is 


a  =  2RXZ. 


(2A) 


(3A) 
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Thus  Che  volume  of  ambient  air  flowing  into  the  channel  per  unit  time  is 


dV 

dt 


T 

mix 


x  2R^Z 


2 

4R  Z 
"t 

mix 


(4A) 


We  now  estimate  che  volume  increase  of  the  channel  (AV  , )  during 

cool 

the  cooling  time  interval  (at  ,)  by  integrating  (4A)  over  At  , .  This 

cooi  '  cool 

gives  us 


AV 


cool 


dV 

dt 


At 

cool 


4R12Z 

"t 

mix 


cool 


(5A) 


The  original  channel  volume  is  given  by 


VQ  =  n  R^Z 


so  that 


AV 


cool 


4  “Ccool 


'cool 


mix 


mix 


(6A) 


(  7a) 


This  relation  is  independent  of  initial  channel  size,  and  the 
volume  increase  will  be  primarily  determined  by  the  channel  temperature 
just  before  mixing  begins  and  Che  ambient  temperature  T^  .  For  both  lightning 
and  the  laboratory  discharges  T^  ~  3000-5000  K  and  T^  ~  300  K.  Thus 
aV  At 

cool  ,  cool  ,  .  . 

— — -  and  — — -  are  about  the  same  in  the  two  cases.  So  our  use  ot 

0  mix 


At  20  T  ,  for  lightning  appears  to  be  justified, 

cool  mix  °  J 

We  note  that  (1A)  can  be  restated  as 


AV 

At 


cool 

cool 


T 

mix 


(BA) 
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The  left  hand  side  is  just  the  average  "mixing  rate"  for  the  problem. 


Thus  the  factor  indeed  appears  to  be  a  fundamental  time  scale  for 

convective  channel  cooling.  Further  investigation  of  the  physical  signif¬ 
icance  of  T  .  should  result  in  a  more  elegant  proof  of  (11) . 
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